Rapid progress on the vertebrate tree of life

<p>Abstract</p> <p>Background</p> <p>Among the greatest challenges for biology in the 21st century is inference of the tree of life. Interest in, and progress toward, this goal has increased dramatically with the growing availability of molecular sequence data. However, we have very little sense, for any major clade, of how much progress has been made in resolving a full tree of life and the scope of work that remains. A series of challenges stand in the way of completing this task but, at the most basic level, progress is limited by data: a limited fraction of the world's biodiversity has been incorporated into a phylogenetic analysis. More troubling is our poor understanding of what fraction of the tree of life is understood and how quickly research is adding to this knowledge. Here we measure the rate of progress on the tree of life for one clade of particular research interest, the vertebrates.</p> <p>Results</p> <p>Using an automated phylogenetic approach, we analyse all available molecular data for a large sample of vertebrate diversity, comprising nearly 12,000 species and 210,000 sequences. Our results indicate that progress has been rapid, increasing polynomially during the age of molecular systematics. It is also skewed, with birds and mammals receiving the most attention and marine organisms accumulating far fewer data and a slower rate of increase in phylogenetic resolution than terrestrial taxa. We analyse the contributors to this phylogenetic progress and make recommendations for future work.</p> <p>Conclusions</p> <p>Our analyses suggest that a large majority of the vertebrate tree of life will: (1) be resolved within the next few decades; (2) identify specific data collection strategies that may help to spur future progress; and (3) identify branches of the vertebrate tree of life in need of increased research effort.</p

Ancestral Vascular Lumen Formation via Basal Cell Surfaces

The cardiovascular system of bilaterians developed from a common ancestor. However, no endothelial cells exist in invertebrates demonstrating that primitive cardiovascular tubes do not require this vertebrate-specific cell type in order to form. This raises the question of how cardiovascular tubes form in invertebrates? Here we discovered that in the invertebrate cephalochordate amphioxus, the basement membranes of endoderm and mesoderm line the lumen of the major vessels, namely aorta and heart. During amphioxus development a laminin-containing extracellular matrix (ECM) was found to fill the space between the basal cell surfaces of endoderm and mesoderm along their anterior-posterior (A-P) axes. Blood cells appear in this ECM-filled tubular space, coincident with the development of a vascular lumen. To get insight into the underlying cellular mechanism, we induced vessels in vitro with a cell polarity similar to the vessels of amphioxus. We show that basal cell surfaces can form a vascular lumen filled with ECM, and that phagocytotic blood cells can clear this luminal ECM to generate a patent vascular lumen. Therefore, our experiments suggest a mechanism of blood vessel formation via basal cell surfaces in amphioxus and possibly in other invertebrates that do not have any endothelial cells. In addition, a comparison between amphioxus and mouse shows that endothelial cells physically separate the basement membranes from the vascular lumen, suggesting that endothelial cells create cardiovascular tubes with a cell polarity of epithelial tubes in vertebrates and mammals

Acoel Flatworms Are Not Platyhelminthes: Evidence from Phylogenomics

Acoel flatworms are small marine worms traditionally considered to belong to the phylum Platyhelminthes. However, molecular phylogenetic analyses suggest that acoels are not members of Platyhelminthes, but are rather extant members of the earliest diverging Bilateria. This result has been called into question, under suspicions of a long branch attraction (LBA) artefact. Here we re-examine this problem through a phylogenomic approach using 68 different protein-coding genes from the acoel Convoluta pulchra and 51 metazoan species belonging to 15 different phyla. We employ a mixture model, named CAT, previously found to overcome LBA artefacts where classical models fail. Our results unequivocally show that acoels are not part of the classically defined Platyhelminthes, making the latter polyphyletic. Moreover, they indicate a deuterostome affinity for acoels, potentially as a sister group to all deuterostomes, to Xenoturbellida, to Ambulacraria, or even to chordates. However, the weak support found for most deuterostome nodes, together with the very fast evolutionary rate of the acoel Convoluta pulchra, call for more data from slowly evolving acoels (or from its sister-group, the Nemertodermatida) to solve this challenging phylogenetic problem

OrthoSelect: a protocol for selecting orthologous groups in phylogenomics

Background: Phylogenetic studies using expressed sequence tags (EST) are becoming a standard approach to answer evolutionary questions. Such studies are usually based on large sets of newly generated, unannotated, and error-prone EST sequences from different species. A first crucial step in EST-based phylogeny reconstruction is to identify groups of orthologous sequences. From these data sets, appropriate target genes are selected, and redundant sequences are eliminated to obtain suitable sequence sets as input data for tree-reconstruction software. Generating such data sets manually can be very time consuming. Thus, software tools are needed that carry out these steps automatically. Results: We developed a flexible and user-friendly software pipeline, running on desktop machines or computer clusters, that constructs data sets for phylogenomic analyses. It automatically searches assembled EST sequences against databases of orthologous groups (OG), assigns ESTs to these predefined OGs, translates the sequences into proteins, eliminates redundant sequences assigned to the same OG, creates multiple sequence alignments of identified orthologous sequences and offers the possibility to further process this alignment in a last step by excluding potentially homoplastic sites and selecting sufficiently conserved parts. Our software pipeline can be used as it is, but it can also be adapted by integrating additional external programs. This makes the pipeline useful for non-bioinformaticians as well as to bioinformatic experts. The software pipeline is especially designed for ESTs, but it can also handle protein sequences. Conclusion: OrthoSelect is a tool that produces orthologous gene alignments from assembled ESTs. Our tests show that OrthoSelect detects orthologs in EST libraries with high accuracy. In the absence of a gold standard for orthology prediction, we compared predictions by OrthoSelect to a manually created and published phylogenomic data set. Our tool was not only able to rebuild the data set with a specificity of 98%, but it detected four percent more orthologous sequences. Furthermore, the results OrthoSelect produces are in absolut agreement with the results of other programs, but our tool offers a significant speedup and additional functionality, e.g. handling of ESTs, computing sequence alignments, and refining them. To our knowledge, there is currently no fully automated and freely available tool for this purpose. Thus, OrthoSelect is a valuable tool for researchers in the field of phylogenomics who deal with large quantities of EST sequences. OrthoSelect is written in Perl and runs on Linux/Mac OS X

The molecular phylogeny of eph receptors and ephrin ligands

<p>Abstract</p> <p>Background</p> <p>The tissue distributions and functions of Eph receptors and their ephrin ligands have been well studied, however less is known about their evolutionary history. We have undertaken a phylogenetic analysis of Eph receptors and ephrins from a number of invertebrate and vertebrate species.</p> <p>Results</p> <p>Our findings indicate that Eph receptors form three major clades: one comprised of non-chordate and cephalochordate Eph receptors, a second comprised of urochordate Eph receptors, and a third comprised of vertebrate Eph receptors. Ephrins, on the other hand, fall into either a clade made up of the non-chordate and cephalochordate ephrins plus the urochordate and vertebrate ephrin-Bs or a clade made up of the urochordate and vertebrate ephrin-As.</p> <p>Conclusion</p> <p>We have concluded that Eph receptors and ephrins diverged into A and B-types at different points in their evolutionary history, such that primitive chordates likely possessed an ancestral ephrin-A and an ancestral ephrin-B, but only a single Eph receptor. Furthermore, ephrin-As appear to have arisen in the common ancestor of urochordates and vertebrates, whereas ephrin-Bs have a more ancient bilaterian origin. Ancestral ephrin-B-like ligands had transmembrane domains; as GPI anchors appear to have arisen or been lost at least 3 times.</p

Acoelomorpha: earliest branching bilaterians or deuterostomes?

The Acoelomorpha is an animal group comprised by nearly 400 species of misleadingly inconspicuous flatworms. Despite this, acoelomorphs have been at the centre of a heated debate about the origin of bilaterian animals for 150 years. The animal tree of life has undergone major changes during the last decades, thanks largely to the advent of molecular data together with the development of more rigorous phylogenetic methods. There is now a relatively robust backbone of the animal tree of life. However, some crucial nodes remain contentious, especially the node defining the root of Bilateria. Some studies situate Acoelomorpha (and Xenoturbellida) as the sister group of all other bilaterians, while other analyses group them within the deuterostomes which instead suggests that the last common bilaterian ancestor directly gave rise to deuterostomes and protostomes. The resolution of this node will have a profound impact on our understanding of animal/bilaterian evolution. In particular, if acoelomorphs are the sister group to Bilateria, it will point to a simple nature for the first bilaterian. Alternatively, if acoelomorphs are deuterostomes, this will imply that they are the result of secondary simplification. Here, we review the state of this question and provide potential ways to solve this long-standing issue. Specifically, we argue for the benefits of (1) obtaining additional genomic data from acoelomorphs, in particular from taxa with slower evolutionary rates; (2) the development of new tools to analyse the data; and (3) the use of metagenomics or metatranscriptomics data. We believe the combination of these three approaches will provide a definitive answer as to the position of the acoelomorphs in the animal tree of life

The enigmatic mitochondrial genome of Rhabdopleura compacta (Pterobranchia) reveals insights into selection of an efficient tRNA system and supports monophyly of Ambulacraria

<p>Abstract</p> <p>Background</p> <p>The Hemichordata comprises solitary-living Enteropneusta and colonial-living Pterobranchia, sharing morphological features with both Chordata and Echinodermata. Despite their key role for understanding deuterostome evolution, hemichordate phylogeny is controversial and only few molecular data are available for phylogenetic analysis. Furthermore, mitochondrial sequences are completely lacking for pterobranchs. Therefore, we determined and analyzed the complete mitochondrial genome of the pterobranch <it>Rhabdopleura compacta </it>to elucidate deuterostome evolution. Thereby, we also gained important insights in mitochondrial tRNA evolution.</p> <p>Results</p> <p>The mitochondrial DNA of <it>Rhabdopleura compacta </it>corresponds in size and gene content to typical mitochondrial genomes of metazoans, but shows the strongest known strand-specific mutational bias in the nucleotide composition among deuterostomes with a very GT-rich main-coding strand. The order of the protein-coding genes in <it>R. compacta </it>is similar to that of the deuterostome ground pattern. However, the protein-coding genes have been highly affected by a strand-specific mutational pressure showing unusual codon frequency and amino acid composition. This composition caused extremely long branches in phylogenetic analyses. The unusual codon frequency points to a selection pressure on the tRNA translation system to codon-anticodon sequences of highest versatility instead of showing adaptations in anticodon sequences to the most frequent codons. Furthermore, an assignment of the codon AGG to Lysine has been detected in the mitochondrial genome of <it>R. compacta</it>, which is otherwise observed only in the mitogenomes of some arthropods. The genomes of these arthropods do not have such a strong strand-specific bias as found in <it>R. compacta </it>but possess an identical mutation in the anticodon sequence of the tRNA_Lys.</p> <p>Conclusion</p> <p>A strong reversed asymmetrical mutational constraint in the mitochondrial genome of <it>Rhabdopleura compacta </it>may have arisen by an inversion of the replication direction and adaptation to this bias in the protein sequences leading to an enigmatic mitochondrial genome. Although, phylogenetic analyses of protein coding sequences are hampered, features of the tRNA system of <it>R. compacta </it>support the monophyly of Ambulacraria. The identical reassignment of AGG to Lysine in two distinct groups may have occurred by convergent evolution in the anticodon sequence of the tRNA_Lys.</p

Evolution of a Core Gene Network for Skeletogenesis in Chordates

The skeleton is one of the most important features for the reconstruction of vertebrate phylogeny but few data are available to understand its molecular origin. In mammals the Runt genes are central regulators of skeletogenesis. Runx2 was shown to be essential for osteoblast differentiation, tooth development, and bone formation. Both Runx2 and Runx3 are essential for chondrocyte maturation. Furthermore, Runx2 directly regulates Indian hedgehog expression, a master coordinator of skeletal development. To clarify the correlation of Runt gene evolution and the emergence of cartilage and bone in vertebrates, we cloned the Runt genes from hagfish as representative of jawless fish (MgRunxA, MgRunxB) and from dogfish as representative of jawed cartilaginous fish (ScRunx1–3). According to our phylogenetic reconstruction the stem species of chordates harboured a single Runt gene and thereafter Runt locus duplications occurred during early vertebrate evolution. All newly isolated Runt genes were expressed in cartilage according to quantitative PCR. In situ hybridisation confirmed high MgRunxA expression in hard cartilage of hagfish. In dogfish ScRunx2 and ScRunx3 were expressed in embryonal cartilage whereas all three Runt genes were detected in teeth and placoid scales. In cephalochordates (lancelets) Runt, Hedgehog and SoxE were strongly expressed in the gill bars and expression of Runt and Hedgehog was found in endo- as well as ectodermal cells. Furthermore we demonstrate that the lancelet Runt protein binds to Runt binding sites in the lancelet Hedgehog promoter and regulates its activity. Together, these results suggest that Runt and Hedgehog were part of a core gene network for cartilage formation, which was already active in the gill bars of the common ancestor of cephalochordates and vertebrates and diversified after Runt duplications had occurred during vertebrate evolution. The similarities in expression patterns of Runt genes support the view that teeth and placoid scales evolved from a homologous developmental module

Site-specific time heterogeneity of the substitution process and its impact on phylogenetic inference

<p>Abstract</p> <p>Background</p> <p>Model violations constitute the major limitation in inferring accurate phylogenies. Characterizing properties of the data that are not being correctly handled by current models is therefore of prime importance. One of the properties of protein evolution is the variation of the relative rate of substitutions across sites and over time, the latter is the phenomenon called heterotachy. Its effect on phylogenetic inference has recently obtained considerable attention, which led to the development of new models of sequence evolution. However, thus far focus has been on the quantitative heterogeneity of the evolutionary process, thereby overlooking more qualitative variations.</p> <p>Results</p> <p>We studied the importance of variation of the site-specific amino-acid substitution process over time and its possible impact on phylogenetic inference. We used the CAT model to define an infinite mixture of substitution processes characterized by equilibrium frequencies over the twenty amino acids, a useful proxy for qualitatively estimating the evolutionary process. Using two large datasets, we show that qualitative changes in site-specific substitution properties over time occurred significantly. To test whether this unaccounted qualitative variation can lead to an erroneous phylogenetic tree, we analyzed a concatenation of mitochondrial proteins in which Cnidaria and Porifera were erroneously grouped. The progressive removal of the sites with the most heterogeneous CAT profiles across clades led to the recovery of the monophyly of Eumetazoa (Cnidaria+Bilateria), suggesting that this heterogeneity can negatively influence phylogenetic inference.</p> <p>Conclusion</p> <p>The time-heterogeneity of the amino-acid replacement process is therefore an important evolutionary aspect that should be incorporated in future models of sequence change.</p

Evolution of the eukaryotic ARP2/3 activators of the WASP family: WASP, WAVE, WASH, and WHAMM, and the proposed new family members WAWH and WAML

<p>Abstract</p> <p>Background</p> <p>WASP family proteins stimulate the actin-nucleating activity of the ARP2/3 complex. They include members of the well-known WASP and WAVE/Scar proteins, and the recently identified WASH and WHAMM proteins. WASP family proteins contain family specific N-terminal domains followed by proline-rich regions and C-terminal VCA domains that harbour the ARP2/3-activating regions.</p> <p>Results</p> <p>To reveal the evolution of ARP2/3 activation by WASP family proteins we performed a "holistic" analysis by manually assembling and annotating all homologs in most of the eukaryotic genomes available. We have identified two new families: the WAML proteins (WASP and MIM like), which combine the membrane-deforming and actin bundling functions of the IMD domains with the ARP2/3-activating VCA regions, and the WAWH protein (WASP without WH1 domain) that have been identified in amoebae, Apusozoa, and the anole lizard. Surprisingly, with one exception we did not identify any alternative splice forms for WASP family proteins, which is in strong contrast to other actin-binding proteins like Ena/VASP, MIM, or NHS proteins that share domains with WASP proteins.</p> <p>Conclusions</p> <p>Our analysis showed that the last common ancestor of the eukaryotes must have contained a homolog of WASP, WAVE, and WASH. Specific families have subsequently been lost in many taxa like the WASPs in plants, algae, Stramenopiles, and Euglenozoa, and the WASH proteins in fungi. The WHAMM proteins are metazoa specific and have most probably been invented by the Eumetazoa. The diversity of WASP family proteins has strongly been increased by many species- and taxon-specific gene duplications and multimerisations. All data is freely accessible via <url>http://www.cymobase.org</url>.</p